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Abstract. Consider n points on the unit 2-sphere. The potential of the interaction of 
two points is a function f(r) of the distance r between the points. The total energy £ of n 
points is the sum of the pairwise energies. The question is how to place the points on the 
sphere to minimize the energy S. For the Coulomb potential f(r) — 1/r, the problem goes 
back to Thomson (1904). The results for n < 5 are well known. The case n = 5 turns out 
to be difficult. For n = 5, Dragnev, Legg, and Townsend [2] give a solution of the problem 
for f(r) = — logr known as Whyte's problem. Schwartz [I] gives a rigorous computer-aided 
solution of Thompson's problem. Hou and Shao give a rigorous computer-aided solution for 
/(r) = — r. We give a solution for biquadratic potentials. 
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1 Introduction 

Let P = {pi, ...,p n } be a configuration of n points on the unit sphere S 2 C M 3 . Let / be 
a real function that we regard as the potential of the interaction of two points. Then the 
potential energy of P has the form 

£f(P) = J2f(\Pi-Pj\), 

i<j 

here \x\ denotes the standard Euclidean norm of x 6 R 3 . The problem of minimizing £/(P) 
is a subject of hundreds of publications (see, e. g., [1] and references there). For the Coulomb 
potential f(r) = 1/r, the problem has a significance in physics and goes back to Thomson 
(1904). 

We are concerned with the energy minimizing problem for small values of n. To eliminate 
trivial or irrelevant potentials, we assume that the function t h-> f(yi) is strongly convex. 
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The results for n < 5 are simple and well-known. They are special cases of a more gen- 
eral result that the regular simplex is a minimizing configuration in S m ~~ l G IR m for every 
dimension m (see, e. g., pQ). 

The case n = 5 already turns out to be difficult. Dragnev, Legg, and Townsend [2] give 
a solution of the problem for f(r) = — logr known as Whyte's problem. Schwartz [4] gives 
a rigorous computer-aided solution of Thompson's problem for n = 5. Hou and Shao [3] 
give a rigorous computer-aided solution for f(r) = —r, for which the problem is well-known 
in discrete geometry. The results of [3] and [I] involve massive calculations that require a 
computer. In all these three cases, a unique minimizer is the so-called triangular bipyramid 
(TBP) that consists of two antipodal points, say the North and South poles, and three points 
on the equator forming an equilateral triangle. However, following pQ and jl] we point out 
that TBP is not universally optimal for either all inverse power potentials / (r) = r~ a , a > 0, 
or potentials of the form f(r) = (4 — r 2 ) k , k G N. For sufficiently large values of a and k, the 
energy takes a smaller value on another configuration, a regular four-pyramid with square 
base (FP). The latter depends on a parameter, say the pyramid's height, whose value is 
chosen to minimize the energy and it depends on /. 

Clearly, a minimizing configuration must be a critical point for the energy Ef. We con- 
jecture that for n = 5 and almost all potentials /, there are only two nontrivial critical 
configurations: TBP and FP. Here we call a configuration trivial if it has repeated points or 
reduces to lower dimension. It would suffice to prove the conjecture for potentials / which 
are polynomials in r 2 . The case of /(r) = — r 2 is degenerate: it is not strongly convex in r 2 . 
In this case the answer is the following: a configuration P is a minimizer exactly if ^2pi = 0, 
see Lemma \2. II below. We consider the first non-degenerate case f(r) = ar 4 — br 2 + c, which 
we call a biquadratic potential. This function is decreasing and strongly convex in r 2 if a > 
and b > 8a. Our main result is the following. 

Theorem 1.1 Let n=5 and let f(r) = ar 4 — br 2 + c, a > 0, b > 8a. Then TBP and FPare 
the only two nontrivial critical configurations for Sf and TBPis a unique global minimizer. 

One motivation for this result is of course to obtain a solution for a non-trivial potential 
without massive computer-aided calculations. We point out that known solutions for sev- 
eral potentials may yield a solution for another potential. Indeed, suppose that for several 
potentials /j, it is known that TBP is a minimizer. Let g = c ifi-> c i > 0' approximate / 
in such a way that g(r) < f(r) for all < r < 2, and g matches with / on all distances r 
between the points of TBP, that is, r = y/2, y/3, and 2. Then TBP is a minimizer for /. For 
instance, if TBP is known to be a minimizer for (4 - r 2 ) k for k = 0,1, 2, 3, 5, then usine; the 
above argument one can show that TBP will be a minimizer for both power potentials r~ a 
and — r a for all < a < 2. Our biquadratic potential covers k = 0, 1, 2, but we admit that 
proving the result for k = 3 and 5 would be hard. 

We prove the main result by solving the system for critical configurations. As a byproduct 
we obtain a result on non-degeneracy of a Cauchy type matrix (Theorem 14.11) that may be of 
interest by itself. We use it to eliminate the bulk of non-critical configurations and reduce the 
analysis of the system to what we call Main Special Case. The latter restricts the problem 
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to configurations having three points on a great circle and two other points symmetric about 
that circle. Surprisingly, even in this simple case the system gives us enough trouble that 
in the end we use symbolic computation software. We point out however that the use of a 
computer is light and can be completely eliminated. 

2 Preliminaries 

We first rewrite the energy in terms of the pairwise inner products. Let x, y G S m ~ l , m G N. 
Let t = x ■ y be the Euclidean inner product, and let r = \x — y\. Then r 2 = 2 — It. For a 
function f(r), we define h{t) = f(s/2 — 2t). If / is decreasing on (0,2], then h is increasing 
on [—1, 1). Then we rewrite the energy £f(P) of a configuration P = {pi, . . . ,p n } in S m ~ l 
in the form 

E h (P) = J2h( Pi - Pj ). (1) 

i<j 

For a biquadratic potential after omitting unnecessary constants we have 

h{t) = (t + a) 2 , a > 1. 
Note that for a > 1 the function h is increasing. We recall the following simple fact. 
Lemma 2.1 Let h{t) = t. Then P = {p\, . . . ,p n } is a minimizer for Eh if and only if 

Proof. Follows from E h {P) = Y.i<jPi 'Pi = \{\ YTi=iPi\ 2 -n).U 

Corollary 2.2 Let n = 5, m = 3. It suffices to show that TBPis a unique minimizer for Eh, 
h(t) = (t + 1) 2 . Then it would follow that TBPis a unique minimizer for Eh, h(t) = (t + a) 2 , 
for all a > 1 . 

Proof. Follows from Lemma [2.11 because for TBP we have J^Pi = 0. □ 

For inverse power potentials f(r) = r~ a , a minimizer can not include repeated points 
because /(0) = oo. The latter is not the case for biquadratic potentials. For convenience we 
consider the case of repeated points separately. 

Proposition 2.3 A minimizer for the biquadratic energy Eh, h(t) = (t + 1) 2 , of 5 points on 
S 2 can not have repeated points. 

Proof. We use the notation E(P) = E h (P) for h(t) = (t + 1) 2 . For TBP the inner products 
are the following: one -1, six 0, and three -1/2. Hence E(TBP) = + 6 + 3(l/2) 2 = 27/4. 

We also use a regular tetrahedron T on the sphere. We recall that T is a universal 
minimizer, in particular, for our potential. For T, all the four inner products equal —1/3. 
Hence E(T) = 6(2/3) 2 = 8/3. 
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Arguing by contradiction, suppose P = {p,p, q\, q 2 , q 3 } is a minimizer. Then we have 

3 

£(TBP) > E{P) = E({p, p}) + E({p, Ql , q 2 , q 3 }) + ^ E({p, %}). 

i=i 

Since E({p,q 1 ,q 2 ,q 3 }) > E(T), then 
3 

E dP> ft» ^ -K(TBP) - £(T) - £({p, p}) = 27/4 - 8/3 - 4 = 1/12. 

i=i 

Hence E({p, $}) < 1/12, z = 1,2,3. Put t = p-q { . Then (t + 1) 2 < 1/12, and t < \/3/6-l < 
—y/2/2. Since the angle between — p and (ft is less that 7r/4, then the angle between (ft and 
5j is less than 7r/2, that is, (ft • > 0, and £7({(ft, > 1. We now have 

27/4 = £(TBP) > E(P) > E({p, p}) + E({q u q 2 , q 3 }) > 4 + 3 = 7, 

which is absurd. □ 

3 Equilibrium equations 

We find the minimizer by finding all critical points of the energy functional. We call them 
critical or equilibrium configurations. 

Let P = {pi, ...,p n } C 5" m_1 . We later restrict to n = 5, m = 3. Put 6^ = pi ■ pj. By 
differentiating the energy (JTJ), using the Lagrange multipliers /Xj, we obtain the equations 

^2 h '{kj)Pj = H&i, l<i<n. 

By multiplying both sides by p,, we can solve for /ij. After eliminating fa we obtain a system 
of equations for critical configurations 

n 

h'ikjXpj - b ijPi ) = 0, 1 < i < n. (2) 

3=1 

Note that if h is good near 1, then we do not have to omit the term with j = i because it 
is equal to 0. The system is difficult to solve in general. Following [2], we point out that in 
the case of the logarithmic potential h(t) = — log(l —t), by adding up all the equations, one 
obtains ^2pi = 0, which is useful in solving the system (J2J). However, for the biquadratic 
potential, this method does not work. Although the equation J^pj = is valid for TBP, it 
is not valid for all critical configurations, in particular for FP. 

Note that the system is invariant under rotations of the sphere. To work with independent 
parameters, one needs a normalization. A natural normalization for m = 3 would consist 
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of fixing one point, say at the North pole, and restricting position of the second point to a 
fixed meridian. However, we use a different normalization. 

Let X be the matrix whose rows are the points pi of the configuration P. Our normal- 
ization requires that the matrix X have orthogonal columns. We show that it is always the 
case in a suitable orthogonal coordinate system in R m . 

Lemma 3.1 Let X be a real n x m matrix. Then there is an orthogonal m x m matrix U 
such that the matrix Y = XU has orthogonal columns, that is, Y*Y is diagonal. 

Proof. We have Y*Y = (XU)*(XU) = U*AU, here A = X*X. The matrix A is a 
symmetric m x m matrix, hence it is diagonalizable in a suitable orthonormal basis. That 
is, there is an orthogonal matrix U such that A = U*AU is diagonal as desired. □ 

Note that in the above proof, the diagonal of A consists of the eigenvalues of A. Also 
note rank(A) = rank(A) = rank(X) < m. Finally, we obtain the following. 

Corollary 3.2 Let X be a realnxm matrix, such that A = X*X is diagonal. Let B = XX* 
and let Xi, X m be the diagonal entries of A. Then the columns of X are eigenvectors of B 
with the eigenvalues Ai, X m , that is BX = XA. The other eigenvalues of B are zeros. 

Let P = {pi,...,p n } C S" 1 ' 1 satisfy the above normalization. Let X = (xik) be the 
matrix of the coordinates of all pi. Let B = = XX* have the eigenvalues Ai,...,A m . 
Then we have 

n n 
i=l j=l 

Here k and I run from 1 to m, b~ki is the Kronecker symbol. Since ba = \pi\ = 1, then 
ti(B) = n. Since tr(XX*) = ti(X*X) = £ A fe , then 

m 

A fc = n. 

fc=i 

We now apply our normalization to the biquadratic energy and the equilibrium equations. 
Denote x = ^2 Pi- Taking into account that 52 ■■ bfj = tr(£> 2 ) = 52 X\, we obtain 

E{P) = J2(bi3 + l) 2 = - + ^ - 2n = 2 E X l + l^l 2 + Y - 2n - ( 3 ) 

i<j M=l k=l 

This expression would attain its minimum if x = and Ai = ... = X m — n/m. However, as 
we will see later (Proposition I5.2p . for m = 3 and n = 5 this situation cannot occur, but it 
can occur for m = 3, n > 5. This observation makes the case n = 5 most interesting. 
For m = 3 and n = 5 we use (E]) to eliminate configurations in lower dimension. 

Proposition 3.3 A minimizer for the biquadratic energy of 5 points on S 2 can not lie in a 
2-plane. 
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Proof. Let P be a critical configuration lying in a plane. Then the plane must pass through 
the origin, otherwise P will not be in equilibrium. By plugging n = 5, m = 2, x = and 
\ 1 = \ 2 = 5/2 in (jSJ we obtain E(P) > 35/4, which is greater than E(TBP) = 27/4. □ 

For the biquadratic potential h(t) = (t + l) 2 the system (j2j) takes the form 

n 

/Xkj + - 6ij»ife) = 0. (4) 

3=1 



Denote 



X = (Xi, ...,X m ) = ^ Pit a i = ^2( b ij + 

i=l j=l 

Then the system (J4]) further reduces to 

(aj - A fc )x iA: = x fc , 1 < z < ra, 1 < fe < m. (5) 
Suppose Xfc 7^ for some 1 < k < m. Then by (J3J) for all i we have 

^ifc = — r • (6) 
a, - A fe 

After summing in i and dividing by Xk, we have 

n 

On the other hand, if Xk = 0, then for every % we have either = Xk or xn~ = 0. 

As we can see, it matters whether or not Xk = 0. Accordingly, we divide the analysis of 
the system into cases depending on how many zeros there are among the numbers xj~- We 
use the notation Case N if there are exactly N zeros. 



4 Case and Cauchy type matrix 

We recall the Cauchy matrix A = (a^), = a . * A with distinct aij and Xk- If A is a square 
matrix, then it is known to be nonsingular, that is, det A ^ 0. We need a similar result here. 
For future references we include a slightly more general result than we need. 

Theorem 4.1 Let 1 < m < n. Let ai, 1 < i < n, and Xk, 1 < k < m, be distinct complex 
numbers satisfying 0. Let A = (a^) be the matrix with entries 



a ik 



t — if k < m, 
a\' m - 1 if k>m. 



Then det A=£0. 
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Proof. Suppose there are numbers such that ^™ =1 r^ik = for all 1 < k < n. We will 
prove all = 0. Introduce 

n n _ 

- A ) 2 tr «* - A 

Then /(Afc) = g(Afc) = for 1 < A; < m. For k > m, we have 

n 

$>of— - l = 0. (8) 

Note that g = reduces to an algebraic equation of degree n, hence it has n roots counting 
multiplicities. The quotient h — f/g is a rational function vanishing at infinity. It has first 
order poles at at and may have poles at the zeros fij of g other than A^. Hence there exist 
numbers Sj and Cj such that 

Here Uj is the multiplicity of fij as a zero of g if //j is not among the numbers A&, and one 
unit less otherwise. Then v j = n ~ m - 

By passing to the limit as A — > ati in the equation (jUJ), we obtain = ?v We will show 
that all Cj = 0. Denote 0(A) = ^ILi ^Za- Let p = n — m. The Laurent expansion of at 
infinity has the form 

-0(A) = A -1 r i + A~ 2 + ... + A-^r.af- 1 + 0(\~ p ~ l ). 

By © we obtain 0(A) = 0(\~ p ~ l ) as A -> oo. Note that 0' = /. Then /(A) = 0(A- p - 2 ), 
and in turn h{\) = 0(A~ p - 2 ). By © we have 

By reducing to common denominator we obtain if) = Q/R, where Q and R are polynomials 
with degi? = Yl Vj = V- Then by (JTOD, Q(A) = 0(A _1 ). Hence ^ = 0, and ft. = 0. 

Since ft/ = /, then the function h satisfies the differential equation h! = gh. Solving this 
equation yields 

h(X) = Cef gWdX = Ce~ 

which is rational only if C = 0. Hence h = and all = as desired. The proof is complete. 

In conclusion we note that the proof simplifies a little for real «j because in this case all 
the zeros of g are simple. The proof further simplifies in the case n — m. □ 

In the case m < n, since the columns of A are linearly independent, and one of them 
consists of units, then we obtain the following. 
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Corollary 4.2 In the assumptions of Lemma 4-l> let m < n. Then there do not exist 
numbers Ck, 1 < k < m, such that Y^k=i ( a - k \ k ) 2 = ^ f or a ^ ^ — * — n ^ 

We now turn to Case 0. 

Proposition 4.3 Let m < n. Then in Case there are no critical configurations without 
repeated points. 

Proof. Let X = (xik) be the matrix of a critical configuration P without repeated points. 
Let Xk 7^ for all 1 < k < m. Then the equation ([6]) holds for all i and k, and ([7]) holds 
for all k. If oti = ctj for some i 7^ j, then by (|6]) we would have pi = pj. Hence the 
numbers a; are distinct. If \k = \i for some k ^ I, then Ylii x ik x ii — Y^i ( a X -\ l k ) 2 ^ 0> so ^ ne 
columns of X would not be orthogonal. Hence the numbers Xk are also distinct. By ([7]), the 
numbers and Xk satisfy the hypotheses of Corollary 14.21 Since \pi\ = 1, then by ([6]) we 

-2 

have J2k (a -Afc) 2 = ^- then Ck = x\ satisfy the equation in Corollary 14. 2\ which is not 
possible. □ 



5 Cases 1, 2, and 3 

We restrict to n = 5, m = 3. We look for nontrivial critical configurations. We regard 
configurations with repeated points and those lying in a plane as trivial. In addition to 
Cases to 3 we introduce Main Special Case, in which three points lie on a great circle, 
call it equator, and the other two are symmetric about the equator. We will come to the 
following conclusion. 

Proposition 5.1 All nontrivial critical configurations in Cases 1, 2, and 3 fall into Main 
Special Case. 

Case 1. Assume for a nontrivial critical configuration, x\ = 0, x 2 7^ 0, and X3 7^ 0. 
Define ii = {i : «j = Ai}. We claim = 2. 

Let i be the complement of I\. First observe that for i G Iq we have Xn = because 
(«j — Ai)xji = Xi = 0. Then I\ 7^ 0, otherwise the configuration lies in the plane X\ = 0, 
hence trivial. Having exactly one nonzero component among xn is not possible either because 
x\ = 0. Hence |ii| 7^ 1. 

Consider |iJ > 3. Since Xk 7^ for k = 2,3, then x-ik = x \ for k = 2,3. Then Xik 
is independent of i G 1% for k = 2,3, so the points pi for i G l\ can differ only in the first 
component Xn. But since \pi\ = 1, then the only remaining freedom is in the sign of xn. So 
if |A| — 3, then the points have to repeat. 

Hence |ii| = 2, the two points Pi, i G I\, are symmetric about the plane X\ = 0, and the 
three points Pi, i G Iq, are in the plane X\ — 0. □ 

Case 2. Let X\ = 0, £2 = 0, and £3 7^ 0. Define I\ and Jo as above. 
Subcase Ai 7^ A 2 . If i G Iq, then xn = 0. If i G ii, then x i2 = because 7^ A 2 . Then 
I ii| = 2 by the same reason as in Case 1, and the configuration falls into Main Special Case. 
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Subcase Ai = A 2 =: A. If i G Iq, then Xn = x i2 = 0, and x i3 = ±1. If i G ii, then 
x i3 = aZa^ i s independent of i. 

If Jo = 0, then all the points pi lie on the same circle, which we regard as trivial. 

If |/o| = 1, then one point is, say (0,0, 1), the North pole, and the rest are on the same 
parallel. Such a configuration will have to be FP, a regular pyramid with square base, which 
falls into Main Special Case because it is symmetric about the plane through the vertex and 
a diagonal of the base. 

If |/o| — 2, then two points are (0,0, ±1), and the rest lie on the same non-equatorial 
parallel. Such a configuration is not critical. 

The case |/ | > 2 can not occur because otherwise the points will repeat. □ 

Case 3. We have x± = x 2 = x 3 = 0. Define Ik = {i : on = A^}, k = 1,2,3. Note that 
every index 1 < i < 5 has to be in one of Ik, otherwise Pi = 0. 

Subcase A&, 1 < k < 3, are distinct. By the same arguments as above we obtain that 
either Ii — or = 2 and the corresponding points pi are (±1,0,0). Similarly, for all k 
the cardinalities \Ik\ can be only or 2, which is not possible because 5 is an odd number. 

Subcase Ai = A2 7^ A3. Then for i G I± = I 2 , we have 2^3 = 0. Also I/3I = 2, and for 
% G J 3 , we have Pi = (0,0, ±1). This configuration will have to be TBP, and it falls into 
Main Special Case. 

Finally, the subcase Ai = A2 = A3 can not occur, which we prove below. □ 

Proposition 5.2 For m = 3, n = 5 there is no configuration for which x = and \± = 
A 2 = A3. 

Proof. Let X be the matrix of the configuration. Then the matrix Y = X/\f\[ has 
orthonormal columns. Add a fourth column to Y with entries all equal to l/y/b. Since 
x = 0, then the resulting matrix still has orthonormal columns. Add the fifth column (y i5 ) 
to obtain an orthogonal matrix. Note the rows of an orthogonal matrix are unit vectors. 
Since all \pi\ = 1, and the lambdas are equal, then |y i5 | will have to be equal for all i, in fact 
?/i5 = ±l/y/E. This column also must be orthogonal to the fourth column, which has equal 
entries, hence = 0- However, the five terms ±l/\/5 can not add up to zero. □ 

Remark 5.3 It turns out that for m = 3, n > 5 the case x = 0, Ai = A 2 = A3 can occur and 
provides minimizers for the biquadratic energy. The minimizers form continuous families. 

6 Main Special Case 

While the a-A-notation proved to be useful in obtaining geometric restrictions on the critical 
configurations, we were unable to use it effectively in solving the equilibrium equations in 
this case. We do it here by brutal force. 
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A configuration P in Main Special Case has the form 

Pi = (cos^,, sin^j, 0), 1 < i < 3, 



p a = (r,0,±Vl - r 2 ), a = 4,5. 

Here 0j, 1 < i < 3, and |r| < 1 are real parameters. For r = it is immediate that the only 
nontrivial critical configuration is TBP. Hence we can assume r ^ 0. We have 

kj = Pi ■ Pj = cos 9 i:j , 9 i: j := 9i-9j, 



b ia = r cos 9i, 



Ho 



2r 2 - 1. 



Instead of using the general equilibrium equations (jlj), we derive them anew in this case. 
The biquadratic energy has the form 

3 

E = E(P)= J2( b ij + 1 ) 2 = (cos% + l) 2 + 2^(rcos^ + l) 2 + 4r 4 . 

i<j<5 i<j<3 i=l 

Then the equations dE/d9i = and dE/dr = yield 

3 

^^(cos 9ij + 1) sin 9ij + 2(r cos 9 { + l)r sin 9i = 0, 1 < i < 3 

3=1 
3 

y^(r cos 9j + 1) cos 9j + 4r 3 = 0. 

3=1 

We rewrite the above equations in terms of the complex variables Zj = e * . 

3 / .2 r 2 



Ef|-^)+ 2 Ef7-7)+ 2 ^ 2 -^ 2 )+4r(, 4 -^) = , l<z<3, (11) 

3=1 Z i / j=l Zi <! 

3 3 

- + zj 2 ) + 2 Y^Zj + + 16r3 + 6r = 0. (12) 



Adding up the three equations (fTTj) . assuming r ^ 0, we obtain 
The latter and ( TT2T) yield 

r E^ +2 E^' = r E ;z 7 2 + 2 E z 7 1 = - 8r3 - 3r - ( 13 ) 
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Introduce the power sums s p = Yl z p P = ±1,±2. We formally rewrite each equation in 
( II ip as a linear combination of powers of z,i with coefficients depending on Zj-s, that is, 

Az 2 + B Zi + Czr 1 + Dzf = 0. (14) 

Then the coefficients have the form 

A = s_ 2 + 2r 2 , B = 2s_! + 4r, C = -2s x - 4r, D = —s 2 — 2r 2 . (15) 

The equation ( 114)) is non-trivial: the case A = B = C = D = can not occur at all, and if 
A = 0, = 1, then the numbers 2, can not be distinct - we leave the details to the reader. 

Let crp, 1 < p < 3, be elementary symmetric functions of zi. Then the three distinct 
numbers Zi satisfy both (JT4J) and z 3 — <j\z 2 + a 2 z — cr 3 = 0. Let z be the fourth root of ( 1T4|) . 
which may coincide with one of Zj. Then we have the identity 

Az A + Bz 3 + Cz + D = A(z - z )(z 3 - a x z 2 + a 2 z - a 3 ). (16) 

The goal of these manipulations is converting the system ( TTTHT2]) to the new variables x — Si, 
y = t = <7 3 . We have 

S! = x, s 2 = x 2 -2yt, s_ 1 = y, s_ 2 = y 2 - 2x/t, o x = x, o 2 = yt, a 3 = t. (17) 

Equating the coefficients of powers of z in ffTB"]) yields 



B = 2y + 4r = -A(z + x), = A(z x + yt), (18) 

C = -2x -4r = -A{z yt + t), D = -x 2 + 2yt - 2r 2 = Atz . (19) 

Eliminating zq = —yt/x yields 

A(-yt + x 2 ) + 2xy + Arx = 0, (20) 

A(-y 2 t 2 + xt) - 2x 2 - Arx = 0, (21) 

Ayt 2 - x 3 + 2xyt - 2r 2 x = 0. (22) 



We modify (I2"0l by multiplying by t and adding ( 12"21 to it. We multiply (T221 by y and add to 
( 12T|) . We plug A = — 2x/t + 2r 2 in (f22l) . Finally, we solve each equation for t and obtain 

^ 3x 2 + 2r 2 x 2 y + 2r 2 y + Ax + 4r ^ 2 x(x 2 + 2r 2 ) . . 

xy 2 + 2r 2 x + 4y + 4r 3y 2 + 2r 2 y(y 2 + 2r 2 ) 

We note that neither numerators nor denominators in ( [23]) can vanish. We already discarded 
r = 0. The case x = can not occur because then y — x — 0, and by (I18H19P we eventually 
get r = ±2, which is not allowed. The other numerators or denominators can not vanish 
either by similar reasons - we leave the details to the reader. 
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Eliminating t yields 

(3x 2 + 2r 2 ) (3y 2 + 2r 2 ) - (x 2 y + 2r 2 y + Ax + Ar) (xy 2 + 2r 2 x + % + 4r) = 0, (24) 
x(x 2 + 2r 2 )(3y 2 + 2r 2 )(xy 2 + 2r 2 x + Ay + 4r) 

- y (y 2 + 2r 2 ) (3x 2 + 2r 2 ) (a; 2 y + 2r 2 y + Ax + Ar) = 0. (25) 

The equations ( 1241125]) reflect only ( II ip . Converting (1T3]) to the new variables yields 

(3x 2 + 2r 2 )(ry 2 + 2y + 8r 3 + 3r) - 2rx(xy 2 + 2r 2 x + Ay + Ar) = 0, (26) 

and the one with the interchange x -<->■ y. The system (124H261) is simple enough for symbolic 
computations software. We consider two cases x ^ y and x = y. 

Let x 7^ y. Note (1231) is antisymmetric in x and y. Then dividing ( 125]) by r(x — y) yields 

8r 4 + Ar 5 x + Ar 2 x 2 + 2r 3 x 3 + Ar 5 y - 8r 2 xy - 8rx 2 y + 2r 3 x 2 y + Ar 2 y 2 

-8rxy 2 + 2r 3 xy 2 + 6x 2 y 2 + rx 3 y 2 + 2r 3 y 3 + rx 2 y 3 = 0. (27) 

Eliminating x and y from (|24| |26| |27|) and (|26|) with the interchange x -H- y using Macaulay2 
yields 

r 4 (l + 3r 2 )(45 + 12r 2 + 2r 4 )(ll + 22r 2 + 12r 4 ) 

x (15645 + 79500r 2 + 159482r 4 + 142416r 6 + 42624r 8 - 10368r 10 + 2592r 12 ) = 0. 

The only real solution is the trivial r = 0, because 42624 — 10368r 2 + 2592r 4 > 0. 

We now consider the case x = y. Recall y = x. Then x = y implies x G R. The 
last equation in (|23|) implies t = ±1. Note that the system in either (zi,r) or (x,y,t,r) is 
invariant under changing the sign of all the variables, so without loss of generality t — 1. 
We have a± = a 2 = x, cr 3 = 1. Then Zj-s are solutions of the equation z 3 — xz 2 + xz — 1 = 
(z - l)(z 2 + (1 -x)z + 1) = 0. Then, say z x = 1, Rez 2 = Rez 3 = (x - l)/2. 

For x = y the equation ( 125]) holds automatically, and the equations ( 124] 126]) for £ = 1 
reduce to 

(3x 2 + 2r 2 ) - (x 3 + 2r 2 x + Ax + Ar) = 0, (28) 
(rx 2 + 2x + 8r 3 + 3r) - 2rx = 0. (29) 

Eliminating x yields 

r(l + 2r)(l + 2r + 2r 2 )(8 - 9r + 6r 2 )(l + 3r + 6r 3 ) = 0. 

The real nonzero solutions are r = —1/2 and the solutions of the equation 

6r 3 + 3r + 1 = 0. (30) 

Let r = —1/2. Then one can find x = 1, Kez 2 = Kez^ = 0, z 2 , z 3 = ±z. This 
configuration is TBP. 
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The equation (I30p has only one real root r —0.286. Eliminating r from ( I281l3"0l yields 
3x 3 — 9x 2 + 15x — 5 = 0. The latter matches with ( 13 Op if r = (x — l)/2. Then the 
configuration consists of (1,0,0), (r, ±vl — r 2 ), 0), and (r, 0, ±yl — r 2 ). This is a four- 
pyramid with square base FP. The inner products have the following values: four r, four r 2 , 
and two 2r 2 - 1. Then the energy is equal to E(FP) = 4(r + l) 2 + 4(r 2 + l) 2 + 2(2r 2 ) 2 w 7.9, 
which is bigger than E(TBP) = 6.75. Hence TBP is a unique minimizer. The proof of 
Theorem 11.11 is complete. 
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